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ABSTRACT 


The  effect  of  in-plane  boundary  conditions  on  the  buckling  loads  of 
simply  supported  ring-stiffed  cylindrical  shells  is  studied,  As  in  the 
case  of  unstiffened  shells,  the  "weak"  in-plane  boundary  conditions  SSI 
and  SS2  yield  here  critical  loads  about  one  half  of  the  "classical" 
loads.  It  was  observed  that  the  SSI  critical  loads  are  identical  with 
the  SS2  loads  and  the  SS4  loads  are  almost  the  same  as  the  "classical" 
SS3  loads. 

The  combined  effect  of  stiffener  parameters  and  in-plane  boundary 
conditions  is  studied.  For  internally  stiffened  shells  the  influence 
of  in-plane  boundary  conditions  is  found  to  diminish  with  increasing 
values  of  stiffener  eccentricity  and  area,  No  such  effect  is  observed 
for  externally  stiffened  shells.  The  buckling  inodes  are  also  studied 
and  found  that  they  are  dependent  upon  shell  length  (or  Z)  and  upon 
stiffener  location  and  parameters, 
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-  coefficients  of  additional  displacements  -  Eqs.  (9) 

-  coefficients  of  displacements  -  Eqs.  (6) 

-  coefficients  of  additional  displacements  -  Eqs.  (13) 

-  coefficients  of  additional  displacements  -  Eqs.  (14) 

-  defined  by  Eqs.  (23) 

-  coefficients  of  axisymmetric  displacement  -  Eqs. (25) 

-  cross  sectional  area  of  ring 

-  distance  between  rings  (see  Fig.  1) 

-  defined  by  Eqs.  (7) 

-  Eh3/12(l-v2) 

-  defined  by  Eqs.  (7) 

-  moduli  of  elasticity 

-  distance  between  centroid  of  stiffener  cross-section 
and  middle  surface  of  shell,  positive  when  inside  (see 
Fig.  1) 

-  shear  moduli 

-  thickness  of  shell 

-  moment  of  inertia  of  ring  cross-section  about  its 
centroidal  axis 

-  moment  of  inertia  of  ring  cross-section  about  the 
middle  surface  of  the  shell. 

-  torsion  constant  of  ring  cross-section 
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-  (1  +  u2^ 

-  length  of  shell  between  bulkheads 

-  CL/2) 

-  moment  resultants  acting  on  element  -  Eqs.  (2) 

-  integer 

-  membrane  forces  resultants  acting  on  element 

-  integer,  also  number  of  half  axial  waves 

-  axial  load 

~  defined  by  Eqs.  (23) 

-  radius  of  shell 

-  defined  by  Eqs.  (23) 

-  defined  by  Eqs.  (23) 

-  defined  by  Eqs.  (23) 

-  number  of  circumferential  waves 


U(n,t)  -  uefined  by  Eqs.  (23) 

u*,  v*,  w*  -  displacements  (see  Fig.  1) 

*  *  * 

u,  v,  w  -  non-dimensional  displacements  (  ■  u  /R;  v  /R;  v  /R 

respectively) 

«0(*)»  vQ(x) ,  wQ(x)  -  additional  displacements  -  Eqs.  (6) 

(uo)n»(v0)n  -  displacements  defined  by  Eqs.  (9) 

x*,  y  *  t  *  *  -  coordinates  (See  Fig.  1) 

*  *  .  * 

x,  y,  z  -  non-dimensional  coordinates  (»  x  /R;  y  /R;  z  /R)  . 

-  (l-v2)1/2(L/R)2(R/h) 


Z 


exponents  -  Eqs.  (9) 

roots  of  characteristic  equation  -  defined  by  Eqs.  (12) 
irR/L 

Kronecker  delta 
(E2A2e2  R/aO) 

(E2Io2/aO) 

(G2It2/aO) 

coefficients  defined  by  Eq.  (14) 

(PR/nD) 

(K3/D)p 

(1-v2)  (E^/E^) 

(1-v2) (E2A2e2/EahR) 


Subscripts  following  a  comma  indicate  differentiation. 
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I .  INTRODUCTION 


The  influence  of  the  in-plane  boundary  conditions  on  the  buckling 
load  of  unstiffened  cylindrical  shells  under  axial  compression  or 
longitudinal  shear  has  been  the  subject  of  many  studies  [1  -  14].  For 
simply  supported  end  conditions,  the  most  significant  result  obtained  was 
that  the  zero  shear  stress  boundary  conditions,  SSI  (Nx  »  N^  =  0)  and 
SS2  (u  «  N  .  «  0) ,  reduce  the  critical  load  to  one  half  of  the  "classical" 
load  -  that  obtained  for  the  SS3  (v  =  =  0)  boundary  conditions,  while 

for  the  SS4  (u  *  v  *  0)  boundary  conditions  critical  loads  equal  to  the 
"Classical"  ones  were  obtained.  For  unstiffened  conical  shells  the  effect 
was  studied  in  [IS],  [16]  and  [17], 


The  influence  of  the  in-plane  end  constraints  on  the  buckling  under 
lateral  and  hydrostatic  pressure  of  unstiffened  cylindrical  shells  was 
investigated  in  [18]  and  [19].  The  effect  in  th6  case  of  free  vibrations 
was  also  studied  in  [20]  and  [21],  The  effect  of  prebuckling  deformations 
on  the  buckling  load  of  unstiffened  shells  was  shown  by  different  investigators 
[9,  20,  23]  to  be  small  except  for  very  short  shells.  For  stiffened  shells 
this  effect  was  found  to  be  even  smaller  [24]  and  [25], 


For  stiffened-cylindrical  shells  the  effect  of  the  in-plane  boundary 
conditions  on  the  buckling  load  was  studied  in  Refs,  [26]  and  [27],  Axi- 
syrometric  buckling  modes  were  not  examined  in  these  studies.  Since,  however, 
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earlier  studies  [28]  showed  that  in  the  case  of  externally  ring-stiffened 
cylindrical  shells  under  axial  compression  these  modes  are  the  critical 
buckling  modes,  they  should  also  be  considered. 

In  the  present  report  a  different  approach  is  employed  to  study  the 
influence  of  the  in-plane  boundary  conditions  on  the  buckling  under  axial 
compression  and  hydrostatic  pressure  of  simply  supported  ring-stiffened 
shells.  The  displacement  method  employed  for  conical  shells  in  [29]  and 
[15]  and  for  cylindrical  shells  in  [14]  (which  is  an  extension  of  that 
used  earlier  for  classical  boundary  conditions  [30]) is  applied.  In  the 
analysis,  a  solution  is  assumed  for  the  w  displacement,  which  solves  the 
first  two  equilibrium  equations  exactly.  The  solution  yields  four  constants 
which  are  determined  by  compliance  with  the  appropriate  boundary  conditions 
(SSi;  2*  3  and  4  )  Then  the  third  equation  is  solved  by  a  standard 
Galerkin  procedure.  Axisymmetric  buckling  modes  are  included  and  an 
extensive  parametric  study  is  carried  out. 

It  say  be  noted  that  the  stiffening  is  assumed  to  b*  closely  spaced 
and  hence  the  stiffeners  are  taken  as  ’’smeared"  or  "distributed"  over  the 
entire  shell,  which  implies  that  discreteness  effects  -  usually  negligible 
[31]  and  [32]  -  are  not  considered.  The  effects  of  eccentricity  of  loading 
are  also  not  included  in  the  present  study. 
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2.  EQUATIONS  AND  BOUNDARY  CONDITIONS 

The  analysis  is  based  on  the  stability  equation  of  [30]  for  buckling 
under  combined  axial  compression  and  hydrostatic  pressure.  For  ring- 
stiffened  cylindrical  shells  these  equations  become: 

Eh/(l-»2)lu>xx.  (i^.)vI4-vw,x]  -  0  U.) 

-(D/R)k2(2*i4|- 

w 

♦  12(R/h)2[(l.u2)(w-v  )  -  vu>x]  ♦  it-jH.)  » 

*  yn2)  *  «.*)>  - 0 

and  che  forces  and  moments  acting  on  an  element  are  given  by: 

Nx  ■  Eh/ (l-v2)  [u  ♦  v(v  ^  -w] 

«  Eh/(l-v2)[(l*y^(v^  -  w)  ♦  vu>x  -  X2*M) 

Um  V  ♦  vfX] 


(lc) 


(2) 


Mx  -  -€»/»)!».„  + 

\  -  -(P/R)Cw>h(1  ♦  no2)  ♦  vw>xx  -  C2(v -  *)3 


M 


(D/R) (1  -  v)w 


”♦*«  .(D/R)[(l  -  V)  ♦  nt2]w  x4 
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In  the  case  of  simply  supported  shells  the  usual  out  of  plane 
boundary  conditions  to  be  satisfied  at  the  edges  are: 

w  ■  M  »  0  at  x  ■  -  (t/R)  or  y  ■  0 

(3) 

x  »  (l/R)  or  y  *  (tt/B) 


The  displacements  u,  v  and  w,  which  are  the  solutions  of  Eqs.  (1),  must 
also  satisfy  the  appropriate  in-plane  boundary  conditions,  one  of  the 
following  4  sets, 


SSI; 

N 

■  N  , 

«  0 

x 

x* 

SS2: 

u 

■  N  , 

*  0 

at  x  *  -(l/R) 

or  y» 

0 

Classical  "SS3; 

V 

«  N 

«  0 

X 

x  «  (t/R) 

or  y  « 

(»/8) 

SS4; 

u 

*  V 

«  0 

«) 

Since  w  *  0  at  the  edges  also  w  *  0  at  the  boundaries  and 
the  condition  for  ■  0  at  the  boundaries  can  be  replaced  by 


w«w*0  at  x  *  -(£/R)  or  y  «  0 

x  -  (t/R)  or  y  «  (*/8) 

On  the  basis  of  [14],  [28]  and  [30]  the  displacements  for  the  solution 


of  Eqs.  (1)  are  chosen  as 


(6) 


The  additional  displacement  wq(x)  is  arbitrarily  assumed  to  be  equal 
to  zero  since  the  Fourier  series  terns  fulfil  the  boundary  conditions  Eqs. 
(S),  Substitution  of  wq(x)  “0  in  the  first  two  stability  equations 
Eqs.  (la)  and  (lb)  yields  two  homogeneous  differential  equations  for  the 
additional  displacements  uQ  and  vq  . 


The  solution  of  Eqs.  (8)  can  be  written  as 


(u^  *  Aeax  sin(t*) 
(v0)n*  **ax  cos(t«) 


(9) 


Substitution  of  Eqs.  (8)  into  (9)  yields  the  following  honogeneous  algebraic 
equation 


r  2 

[o  - 

(■^~)t2]  -  (ij^)at 

r 

A 

e 

f 

_ 1 

((^)a2  -  (1  *  u2)t2 
• 

X  4 

B 

>  - 

and  non  vanishing  values  for  A  and  B  are  obtained  if  the  determinant  of 
the  coefficients  of  A  and  B  is  vanishing.  Thus  the  characteristic 
equation  for  a  is  obtained 


-  t^Kl*^)  -  v]  o2  ♦  t*(l  *  ^(“J”)  *  0 


(ID 


and  the  roots  of  this  equation  are 


w(k  -  v)  ♦  J[k  -  1)  (k  -  v  ), 

alm  ’  a2  *  t[.V - '  S - “I 


1/2 


, 


•3  •  -  «4  -  t\*  - 


(12) 
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wheve  k  «  (1  ♦ 

2 

It  can  be  shown  that  Oj  and  o2  are  always  real,  because  (k-l)(k-v  )>  0. 
Similarly  02  and  are  always  real  for  all  practical  applications  because 

(k  -  v)  >|{k-  1)  (k  -  v2)  for  a  large  range  of  k. 


With  the  above  calculated  values  of  a,  the  additional  displacements  uQ 
and  v  become 


(u  )  «  C  sin(t$)  I  A,  l 
0  n  n  j.j  Jn 


o^x 


(13) 


<V»"  Cn  co,tt*>  “jn1 


V 


From  Eqs.  (10)  the  relation  between  A^n  and  B^n  is  given  by 


where 


and 


V  ■  V'jn 


2a  -  Cl  -  v)t 
ej  *  U  +  v)o.t 


1 


ei  *  ‘  02  ;  e3  *  ’  04 


(14) 


Hence  the  complete  displacements  can  be  written  as 

m 

u  «  sin(t$)  I  Cn[ancos(n8y)+Alnsh(a1x)+A2nsh(oJx)+A3nch(o1x)+A4nch(o3x)] 


v  «  cos(tf)  E  Cn(bnsin(nSy)+e1Alnch(a1x)+e3A2nch(a3x)+e1A3nsh(o1x)+ 
n»l  1 

♦  e3A4nsK  (ojX)  ] 


Si.;(t4)  Z  C  sin(nBy) 
n-1  n 


w  * 


(15) 
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These  equations  include  four  constants  of  integration  A.  which  will 
be  determined  from  the  appropriate  in-plane  boundary  conditions. 
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3.  COMPLIANCE  WITH  IN-PLANE  BOUNDARY  CONDITIONS 


The  values  of  the  constants  A.  axe  determined  by  enforcing  the  four 

jn 


sets  of  boundary  conditions  listed  in  (4). 


For  example,  in  case  SS4  the  requirement 


u  ■  v  *  0  at  x  »*(t/R)  or  y  *  0 

x  a  (t/R)  or  y  *  (w/B) 


yields 


sh(ctji/R) 

shfajl/R) 

8jCh(aj£/R) 


shCa3£/R) 


sh(Oj£/R) 


-chCOjt/R) 

ch(art/R) 


03ch(a3t/R)  -ejShCajt/R) 


fijChCajl/R)  63ch(a3l/R)  BjShCojt/R) 


-  ch<»3t/R>  1  fAJ  r  *„ 


chfaji/R)  A 


(-l)V 


•63sh(a3l/R)  AJn|  0  j 

03sh(a3i/R)j  [A,  j  [o  j 


This  matrix  equation  can  be  divided  into  two  matrix  equations  one 
yielding  symmetric  modes  of  buckling  -  n  *1,  3,5...  and  one  antisymmetric 
modes  -  n  «  2,4,6... 


For  symmetric  modes: 


n  *  1,  3,  5... 


AJn  ’  *  0 


sh(ojt/R)  shCctjl/R) 


8jCh(oji/R)  83ch(a3l/R) 


(17a) 


-  10  - 


For  antisymmetric  modes: 


*  2,  4 1  6«  *  • 


Aln  '  A2n  ’  0 


chCajt/R)  ch(o3t/R)  A3n 


OjShCo^/R)  63sh(a3t/R) 


(17b) 


In  case  SS3 


N  «  v  «  0 
x 


In  a  similar  manner  to  case  SS4  ,  a  set  of  4  homogereous  equations  is 
obtained  for  which  only  the  trivial  vanishing  solution  exists 

Aln  -  A2n  *  A3n  *  S  '  0  (W) 

and  the  "classical"  solution  of  [28]  is  obtained. 


In  case  SS  2 


u  ■  N  «  0 

For  symmetric  modes:  n  *  1,  5»  5 


s  ■  s  ■ 0 


sh(a,i/R) 


sh(a3l/R) 


OjO^shCa^t/R)  a3e3sh(a3t/R) 


r 

an 

nBbn 


(19a) 


(20b) 
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4.  SOLUTION 


Now  every  term  of  the  displacements  series  Eqs.  (15)  is  a  solution 
of  the  first  two  stability  equations  (la)  and  (lb) . 

The  third  equation  (lc)  is  solved  by  a  Galerkin  procedure. 


2s  U/R);(ir/B) 

fO  d\-  t/R, ;0  -(D/R)U2(2“.H-V,W)*",XXXX*'2*''tl5",XX*4*(1*n02)".W* 


W 


d2(R/h)‘[(l+u2)(w-V^  v)U<x]+X(^H.)+Xp[(-^-+  w^]} 


.XX 


.  wmdx  »  0 


(21) 


Substitution  of  Eqs.  (15)  and  performing  the  necessary  integration 
yields  a  system  of  linear  homogeneous  algebraic  equations  for  the  un¬ 
knowns  X  and  X  which  are  the  critical  load  parameters 
P 


£  C  [A(n,m)  ♦  XB(n,m)  D(n,m)J  =0  (22) 

n 

m  «  1,2,3, , , . 

where  W*  n  *  1,3,  5...  for  symmetric  modes  of  buckling  and 

n  ■  2,4,  6...  for  asymmetric  modes  of  buckling 
and  where  for  the  symmetric  modes 

A(r,,m)  -  6m  ^Q(n,t)  ♦  ch(a1t./R)R(m)AJn  +  ch(a3*/R)S(m)A2n 

(23a) 

B(n,m)  -  6  T  (n)  (23b) 

Hal 

D(n,m)  *  «  U  (n,  t)  (23c? 

antisymmetric 

and  for  the  amjpmmammAe  modes 


A(n,ra) 

8(n,ra) 

D(n,m) 

Q(n,t) 

Q(n.t) 

R(»)  » 


s  ^tt^B  Q(n» -  sh^i/R^fm)  AJr  -  sh (»3 i/R ) S (ra)  A4r 

-  6«,T<»> 

*  <«u(n»t) 

,  R(m),  S(m),  T(n),  U(n,t)  in  Eqs.  (23)  are  defined  as  follows: 

•  {U2t2-12(R/h)2(l*u2)}tbR  -  12(R/h)2tvn8an  +  (l+u2)]-nV 
-  t2(t2(l  ♦  no2)  ♦  (2  ♦  nt2)n2e2  -  2c2J) 

20  UigTt  m  ^  Q* 

"j^—n  +  w2)  -v0n 


(23d) 

(23e) 

(23f) 


20_m8n  -  *  a, 

S(m)  -  {;  t*  -  12(R/h;  [t(l  ♦  M„)  -  vfe=)]> 

ct2  -mb2  8*  2  *  05 

?  2  a2 

T(n)  »  (*72$) 

U(n,t)  *  (ir2/20)[t2  +  S^!.) 

Truncating  of  the  displacement  series  at  n  •  N  yields  a  NxN  stability 
matrix  Eq,  (22)  whose  lowest  eigenvalue  yields  the  critical  load  for  n  »  M. 
in  the  calculations,  the  circumferential  number  of  waves  (t)  is  treated  as  .t 
parameter  for  which  the  minimal  critical  value  is  found  for  a  given  n.  The 
size  NxN  of  the  matrix  is  determined  by  the  convergence  criterion  for  the 
critical  load. 
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5.  AXISYMMETRIC  BUCKLING. 


For  this  mode  of  buckling  t  *  0  and  there  is  no  dependence  on  v 
and  e  and  the  equilibrium  equation  (1)  degenerate  to 


(Eh/l-v2)  vu  -  vw  ]=0 

f  AA  |  A 

■(D/R)['W,xxxx+  12(R/h)2Cvu  x  -  w)  -  (*♦*  H^r)]  =  0 


(24a) 


(24b) 


If  as  before,  wQ  »  0,  the  additional  displacement  uQ  becomes 


"o  ■  'V>Ax  X  *  <\2hx 


and  the  complete  displacements  are  then 
00 

“  ‘  nil  C"(l"  C05l"B!')  *  (A’>1)^*!Ao2)Ax1 


w  *  l  C  sin(nBy) 
n«l  n 

and,  as  before,  the  coefficients  are  determined  by  compliance  with  the 
appropriate  boundary  conditions. 

For  axisymmetric  buckling  the  in-plane  boundary  conditions  (4)  are  given 


In  cases  S.S.l  and  S.S.3:  ■  0  at  x  *  -(l/R)  or  y  *  0 


and  cases  S.S.2  and  S.S.4:  u  ■  0 


x  ■  (t/R)  or  y  «  (ir/B) 


The  coefficients  are  therefore: 
For  6  S  4  and  S.S.2 


(Art1) 


ol'Ax 


—  S 
n  n 


n  ®  1,  3,  5. . . 


n  *  2,  4,  6. . . 


15  - 


(AJ,  =  { 
v  02'' Ax 

-  a 

l  n 


n  =  1.3,5. . . 

^  "  2,4»6i § . 


(28) 


and  for  SS5  and  SSI 

(A0>5AX  "  ^oz’ax’0  (29) 

(Ao2^Ax  is  assumed  to  vanish  because  it  represents  a  rigid  body 
translation  in  the  axial  direction  and  thus  can  be  ignored: 

The  second  stability  equation  (24b)  is  again  solved  by  the  Galerkin 
procedure: 


2"  2  w  M 

/  d^  f  (D/R)  [w  ♦  12(R/h)  (1+p-)  (vu  -w)+- 4— — (A+A  )Jw  dx  *  0 
o  -1/R‘  o  .xxxx  i  ,x  2  P  m 

'  (30) 

which  yields  for  the  symmetric  modes  n  «  1,3,5...  the  set  of  linear  algebraic 
homogeneous  equation 

cn{t*2/23)6m(24(R/h)2[an  nBve(l+v2)]+2  nV]  >+(48/m^CAol)Ax(R/h)2,rv  ♦ 

♦  nV(n2/2B)  S^X  ♦  X) }  *  0  a  »  1,3,5,....  (31) 

and  for  the  asymmetric  modes  -  n  «  2,  4,6..  another  set 

6mnCn{l,2/2e){24(R/h)2fan  n6v  +  (1+y2^  +  ■  +n2e2(X+>.p) }  *  0 


m  *  2,4.6.  ... . 


(32) 
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6.  NUMERICAL  RESULTS  AND  DISCUSSION 

In  the  numerical  work  two  main  stiffener  configurations,  differing 
in  their  eccentricity,  were  studied,  one  with  (e2/h)  =  *1,  (A2/ah)  *  0.5, 

3 

(I22/ah  )  »  2  and  nt2  *  3  (Table  la)  and  the  other  one  with  (e2/h)  =  *5, 

(A2/ah)  =  0.5,  (I20/ah^)  =  2  and  nt2  a  3  (Table  lb).  In  both  studies 
the  shell  geometry  was  varied  in  the  ranges  0.03  ^  (L/R)  ^  2.00  and 
100  <  R/h  ^  2000.  The  critical  loads  obtained  for  these  geometries  of 
shells  and  stiffeners  are  presented  in  Tables  la  and  lb.  These  loads  are 
computed  for  the  four  simple  support  in-plane  boundary  conditions  SSI  to 
SS4.  The  calculated  critical  loads  are  compared  with  the  "classical” 
bucklings  loads  (SS3  boundary  conditions)  and  the  results  are  plotted  as 
a  function  of  the  Batdorf  parameter  Z  in  Figs.  2a  and  2b  and  as  a 
function  of  the  nondimensionalized  shell  length  (L/R)  in  Fig.  3  with  (R/h) 
as  an  additional  parameter. 

From  Figs.  2a  and  2b,  as  well  as  Fig.  3,  it  can  be  seen  that  the  SSI 
boundary  conditions  yield  exactly  the  same  critical  loads  as  the  SS2 
boundary  conditions  and  that  the  SS4  boundary  conditions  yield  the  same 
loads  as  the  "classical"  SS3  conditions.  From  Figs.  2a  and  2b  it  can  be 
seen  that  in  the  very  low  range  of  the  Batdorf  parameter  Z  <  0.2  the 
different  in-plane  boundary  conditions  do  not  differ  the  critical  loads. 

But  from  Tables  la  and  lb  one  may  also  observe  that,  though  the  critical 
loads  are  identical  the  corresponding  buckling  modes  are  completely  different. 
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The  SSI  and  SS2  boundary  conditions  yield  asymmetric  buckling  nodes  where¬ 
as  the  SS3  and  SS4  boundary  conditions  yield  axisymmetric  buckling  modes. 

It  may  be  noted  that  for  asymmetric  modes  the  critical  number  of 
circumferential  waves  is  t  «  2  in  many  cases.  For  such  small  numbers  of 
circumferential  waves  the  Donnell  type  stability  equations  employed  in  the 
analysis  are  inaccurate  in  general.  However,  since  the  critical  loads  are 
found  to  increase  only  very  slightly  if  t  is  increased  from  2  to  4  or  S 
(usually  by  less  than  1%),  the  buckling  loads  with  t  »  2  can  here  be  taken 
as  close  approximations. 

With  increasing  values  of  Z  up  to  values  of  Z  v  4.7  the  critical  loads 
corresponding  to  SSI  and  SS2  boundary  conditions  decrease  to  values  as  low 
as  35%  of  the  "classical"  critical  loads  for  Z23  4.7.  Beyond  this  value 
of  Z,  the  SSI  and  SS2  critical  loads  increase  with  Z,  reaching  values  of 
about  one  half  of  the  "Classical"  load  for  Z  ^  10.  For  Z  >  10  there  is  a 
slight  increase  of  the  SSI  and  SS2  critical  loads  with  increasing  Z.  This 
behavior  differs  slightly  from  that  of  unstiffened  shells  for  which  the  SSI 
and  SS2  boundary  conditions  yield  half  the  critical  "classical"  load  in¬ 
dependent  of  _  (See  [6],  [7],  [13]  and  [14]). 

Figures  2a  and  2b  show  clearly  that  the  influence  of  the  in-plane 
boundary  conditions  depends  on  the  eccentricity  of  the  stiffeners.  In  Fig. 

2a  it  is  observed  that  in  the  case  of  low  values  of  eccentricity,  (ej/h)  «  *1 
both  externally  or  internally  stiffened  shells  are  equally  affected  by  the 
in-plane  boundary  conditions  for  values  of  Z  larger  than  10.  A  noticeable 
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difference  between  external  and  internal  stiffening  appears  there  only  in 
the  range  4  <  Z  <  10  where  first  external  stiffening  yields  lower  ratios  of 
critical  loads  than  internal  stiffening  and  then  the  effect  is  inverted 
but  of  smaller  magnitude. 

With  the  larger  vclues  of  eccentricity  (e2/h)  *  *S  in  Fig.  2b  a  noticeable 
difference  is  evident  between  internal  and  external  stiffening  for  low  value  of  Z 
(Z  >  1)  as  well  as  for  larger  values  of  Z,  External  stiffening  yields  lower 
ratios  of  SSI  and  SS2  to  SS3  loads  for  most  shell  geometries,  except  for  a 
small  region  6  <  Z  <  9  where  the  effect  is  inverted.  Similar  conclusions  can 
be  drawn  from  Fig.  3. 

There  is  some  scatter  of  the  points  representing  the  ratios  of  the 
critical  loads  when  plotted  versus  the  shell  geometry  parameter  Z.  There 
are  two  reason*  for  this  scatter.  First  the  shell  geometry  parameter  is 
only  an  approximate  overall  parameter  in  the  case  of  stiffened  shells  and 
hence  most  of  the  scatter  disappears  when  the  ratios  are  plotted  for  the 
separate  geometric  parameters  (See  Fig.  3) .  Secondly  the  critical  values 
occur  at  different  circumferential  curve  numbers  for  different  in-plane 
boundary  conditions.  The  necessity  of  integral  values  causes  the  well 
known  "ripples"  in  the  curves  of  buckling  load  versus  geometry  parameters, 
which  differ  for  each  case  ai.d  result  in  scatter  when  divided  to  give  the 
ratios  of  the  buckling  loads. 
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In  [28]  it  was  concluded  that  externally  ring-stiffened  cylindrical 
shells  should  always  buckle  in  an  axisymmetric  mode.  The  present  results 
(Tables  la  and  lb)  extend  this  conclusion  in  the  SS  4  boundary  conditions 
which  are  also  shown  to  yield  axisymmetric  modes  of  buckling  for  such  shells. 

In  Table  2  and  Fig. 4  the  variation  of  the  influence  of  the  in-plane 
boundary  conditions  with  increasing  eccentricity  is  studied.  The  geometries 
of  the  shells  are  given  in  Table  2.  From  Fig.  4  it  can  be  seen  that  for 
external  stiffening  the  influence  of  in  plane-boundary  conditions  does  not 
change  with  magnitude  of  the  eccentricity,  for  inside  stiffening,  on  the  other 
hand,  the  magnitude  of  the  eccentricity  affects  the  influence  of  the  in-plane- 
boundary  conditions  which  is  reduced  with  increasing  eccentricity. 

The  variation  of  the  influence  of  the  in-plane  boundary  conditions  with 
increase  of  the  rings  area  parameters  (A^/ah)  is  investigated  in  Table  3 
and  Fig.  5.  Geometry  and  dimensions  of  the  shells  studied  are  presented  in 
Table  3.  From  Fig.  5  it  is  seen  that  for  externally  stiffened  shells  the 
magnitude  of  the  ring-area  is  practically  inmaterial  whereas  for  inside  stiffening 
it  is  a  major  factor.  The  SSI  and  SS2  critical  loads  increase  noticeably 
with  increasing  values  of  the  area  parameter  (A2/ah)  and  approach  the 
"classical" (SS3;  critical  load. 

In  Table  4,  the  effect  of  increasing  the  moment  of  inertia  parameters 
(I22/ah3),  in  relation  to  the  influence  of  the  in-plane  boundary  conditions 
is  studied  on  some  of  the  shells.  The  dimensions  and  geometry  of  the  shells 
examined  are  given  in  Table  4.  The  moment  of  inertia  is  found  to  have  no 
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effect  on  the  critical  load  ratios  of  the  "weak"  boundary  conditions  SSI 
and  SS2.  In  Table  5  and  Fig.  6  the  effect  of  increase  of  moment,  of  inertia 
is  also  studied  for  a  wider  range  of  the  parametei,  1  <  I22/ah5  <  50,  with 
the  same  negative  result.  The  dimensions  and  geometry  of  the  shells  are 
listed  in  Table  5. 

The  buckling  modes  of  the  w  displacement  (radially  inwards)  were 
also  studied  (Tables  la  and  lb).  Some  of  the  results  are  given  in  Figs. 

7  to  11. 

In  Figs.  7a  to  7d  the  buckling  modes  for  "thick"  shells  with  (R/h)=100 
and  low  eccentricity  (e2/'h)  =  *1  are  presented  ,  It  can  be  seen  from  these 
figures  that  short  shells  (Fig.  7a  -  Z  =  .954  and  Fig.  7b  Z  *  8.59  )  yield 
identical  modes  for  all  the  in-plane  of  boundary  conditions  (SSI  to  SS4) . 

Note  also  that  the  modes  are  identical  for  internal  and  external  stiffening. 

As  the  length  of  shells  increases  (Fig.  7c  -  Z  *  95.4  and  Fig.  7d-  Z  =382  ) t 
the  buckling  modes  for  the  SSI  and  SS2  boundary  conditions  remain  identical 
and  independent  or  stiffener  location.  It  can  be  seen  that  these  "weak" 
boundary  conditions  are  characterized  by  edge  buckling  which  is  the  reason 
for  their  low  critical  loads.  On  the  other  hand  the  SS3  and  SS4  boundary 
corditions  yield  modes  that  differ  completely  from  those  of  SSI  and  SS2  and 
depend  also  slightly  on  the  stiffener  location.  Note  that  for  internal 
stiffening  the  SS4  mode  differs  slightly  from  the  SSa  mode  while  for  external 
stiffening  the  modes  are  identical.  Fig.  7d  shows  a  more  pronounced  difference 
between  the  SS3  and  SS4  modes  with  further  increase  of  the  shell  length.  The 
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SS3  mode  now  significantly  differs  from  the  SS4  mode,  and  internal  and 
external  stiffeners  yield  different  modes  of  buckling. 

In  Figs.  8a  to  8d  the  buckling  modes  are  given  for  similar  shells 
but  with  a  larger  ring-eccentricity-  (e2/h)  *  1  5.  Increase  in  eccentricity 
changes  the  buckling  modes.  The  SS3  and  SS4  buckling  modes  are  now  not 
identical  with  the  SSI  and  SS2  modes,  even  for  short  shells:  Fig.  8a-  Z  *  .954 
and  Fig.  8b  -  Z  *  8. 59  ,  For  the  short  shell  (Fig.  8a)  no  dependence 

upon  the  location  of  stiffeners  is  observed  for  all  in-plane  boundary 
conditions,  whereas  location  dependence  can  be  seen  for  the  medium  shell 
of  Fig.  8b.  For  this  shell,  the  SS3  mode  differs  for  internal  and  external 
stiffeners.  With  increase  of  the  shell  length,  the  modes  of  Figs.  8c  and 
8d  are  obtained  which  exhibit  a  behavior  similar  to  that  in  Figs.  7c  and 
7d. 

*  4 

In  Figs.  9a  to  9d  the  buckling  modes  are  presented  for  thin  shells, 

(R/h)  *  2000,  having  a  low  value  of  eccentricity.  The  reduction  of  shell 

thickness,  or  rather  the  inrease  in  Z,  yields  different  modes  for  the  SSI 
and  SS2  boundary  conditions  than  for  SS3  and  SS4.  Note  that  for  the  short 
shell  of  Fig.  9a  -  Z  *  19.1  the  modes  are  location  independent.  With 
increase  in  length  of  shell  (Fig.  9b  -  Z  «  477  )the  SS3  inodes  begin  to 
differ  from  the  SS4  modes,  and  the  SS4  modes  are  also  dependent  upon  the 
location  of  stiffeners.  From  Figs.  9c  and  9d  it  can  be  seen  that  with 
further  increase  in  shell  length  the  SSI  and  SS2  boundary  conditions  yield 
identical  buckling  modes  that  are  location  independent.  Similar  conclusions 
apply  to  the  SS3  and  SS4  boundary  conditions. 
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In  Figs,  10a  to  lOd  the  buckling  inodes  are  studied  for  shells  with 
the  same  geometry  as  Figs. 9a  to  9d,  except  for  larger  eccentricity, 

(®2/h)  *  *5.  For  short  shells,  the  results  of  Fig.  10a  -  Z  =  19.1  and 
Fig.  10b  -  Z  *  477  are  similar  to  those  of  Fig.  9a  and  9b.  With  a 
further  increase  of  shell  length,  a  significant  effect  of  the  increased 
eccentricity  is  noted  for  SS4  modes  of  internally  stiffened  shells  (Figs. 

10c  -  Z  «  1910  and  lOd  -  Z  =  7630).  For  corresponding  externally  stiffened 
shells  no  similar  effect  is  observed. 

The  influence  of  stiffener  eccentricity  on  the  SS4  modes  in  the  case  of 
long  shells  is  studied  in  Figs.  11a  to  lie,  A  significant  effect  is  only 
observed  for  large  values  of  eccentricity  (Figs,  lib  -  1  *  86C0 ;  «2/h  *  5 
and  11c  -z  ■  8600; (e2/h)  •  5).  From  Figs,  lib  and  11c  it  can  also  be  seen 
that  an  increase  in  the  moment  of  inertia  of  the  stiffeners  changes  the 
SS.4  modes  noticeably. 

It  should  be  pointed  out  that  in  all  the  "long"  shells  the  SSI  and  SS2 
modes  are  always  chnracterized  by  edge  buckling,  which  can  explain  the 
significant  reduction  in  the  critical  loads  corresponding  to  these  boundary 
conditions.  These  boundary  conditions  yield  buckling  modes  which  were  almost 
always  one  sided  -  inwards  (positive  w  displacement) >  except  for  long  shells 
with  large  eccentricities  as  in  Figs.  Ub  and  11c. 
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7.  CONCLUSIONS 


a.  As  for  unstiffened  shells,  the  critical  loads  of  ring-stiffened 
shells  depend  on  the  in-plane  boundary  conditions.  The  "weak" 

SSI  and  SS?  condition  yield  identical  buckling  loads  which  are 
about  one  half  the  "classical"  SS3  loads.  The  SS4  boundary  con¬ 
ditions  yield  critical  loads  which  are  practically  equal  to  the 
"classical"  loads,  or  very  slightly  larger. 

b.  The  buckling  loads  of  very  short  shells  -  Z  <  0.1  are  independent 
of  the  in-plane  boundary  conditions. 

c.  In  the  range  of  0.1  <  Z  <  10  the  SSI  and  SS2  boundary  conditions 
yield  critical  loads  which  are  as  low  as  35%  of  the  corresponding 
"classical"  load. 

d.  For  externally  stiffened  shells  the  influence  of  in-plane 
boundary  conditions  is  not  affected  by  the  stiffener  geometry. 

e.  For  internally  stiffened  shell,  on  the  other  hand,  critical  loads 
for  the  "weak"  in-plane  boundary  conditions  SSI  and  SS2  increase 
with  stiffener  area  and  eccentricity.  For  very  large  stiffener 
area  they  approach  the  classical  SS3  critical  loads.  Changes  in 
the  moment  of  inertia  of  stiffeners  have  a  negligible  effect  on  the 
influence  of  in-plane  boundary  conditions. 
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f.  The  buckling  modes  depend  on  the  shell  length  (or  Z)  and  on  the 
stiffeners  geometry. 
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FIG.  5  EFFECT  OF  STIFFENER  AREA  ON  THE  RATIO  OF  BUCKLING  LOADS 
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1  J  ABSTRACT 

The  effect  of  in-plane  boundary  conditions  on  the  buckling  loads  of  simply 
supported  ring-stiffened  cylindrical  shells  is  studied.  As  in  the  erse  of 
unstiffened  shells ,  the  "weak”  in-plane  boundary  conditions  SSI  and  SS2  yield 
here  critical  loads  about  one  half  of  the  "classical"  loads.  It  was  observed 
that  the  SSI  critical  loads  are  identical  with  the  SS2  loads  and  the  SS4  loads 
are  almost  the  same  as  the  "classical"  SS3  loads.  The  combined  effect  of 
stiffener  parameters  and  in-plane  boundary  conditions  is  studied.  For  internally 
stiffened  shells  the  influence  of  in-plane  boundary  conditions  is  found  to 
diminish  with  increasing  values  of  stiffener  eccentricity  and  area.  No  such 
effect  is  observed  for  externally  stiffened  shells.  The  buckling  modes  are  also 
studied  and  found  that  they  are  dependent  upon  shell  length  (or  Z)  and  upon 
stiffener  location  and  parameters. 
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